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We consider an expanding boost-invariant plasma at strong coupling using the AdS/CFT corre- 
spondence for TV = 4 SYM. We determine the relaxation time in second order viscous hydrodynamics 
and find that it is around thirty times shorter than weak coupling expectations. We find that the 
nonsingularity of the dual geometry in the string frame necessitates turning on the dilaton which 
leads to a nonvanishing expectation value for trF 2 behaving like r~ 10 / 3 . 
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I. INTRODUCTION 

It is believed that the quark-gluon plasma (QGP) produced at RHIC in heavy-ion collisions is strongly coupled (see 
e.g. [H) and is well described by almost perfect fluid hydrodynamics [2]. Therefore it is very interesting to develop 
methods for studying its properties nonperturbatively from first principles. A very powerful tool in the study of the 
dynamics of gauge theory at strong coupling is the AdS/CFT correspondence Although so far there does not exist 
a version of the correspondence with a gauge theory which would have all the features of QCD, even the simplest 
version for J\f = 4 Super Yang-Mills theory was argued to share a lot of properties at finite temperature with QCD 
CN| plasma. But of course one has to keep in mind the differences which may not be important for some features of the 
dynamics but which may be crucial for other phenomenae. In this paper we thus work exclusively with hot expanding 
plasma in the Af = 4 SYM theory. 

Extensive work has been done in the understanding of transport properties of the plasma at fixed temperature like 
calculating the shear viscosity [3, H[ ■ Much less is known about the properties of more dynamical time-dependent 
processes. On a qualitative level, thermalization has been suggested to correspond to black hole formation in the bulk 
£NJ , of the five-dimensional dual geometry Q , while cooling was advocated to correspond to black hole motion in the 5th 
■ direction 0. 

In Q a quantitative framework has been proposed for studying the time-dependent expansion of a boost-invariant 
plasma system. The criterion of nonsingularity of the dual geometry was shown to predict almost perfect fluid 
hydrodynamic expansion Q with leading deviations coming from shear viscosity Q with the shear viscosity coefficient 
being exactly equal to the one derived in the static case in [4j. Further work in this framework include Hamuli, 
[l~3l . [3]. The aim of this paper is to investigate in more detail the hydrodynamic expansion and to determine the 
remaining parameter in second order viscous hydrodynamics [lol [l6| - the relaxation time Trj. This requires going 
^ ' one order higher in the subasymptotic expansion of the geometry. 

The plan of this paper is as follows. In section II we will describe the kinematic regime of longitudinal boost 
invariance. Then in section III we will briefly review second order viscous hydrodynamics. In section IV we describe 
the AdS/CFT methods used here and review, in the following section, the results obtained so far and a method of 
determining the relaxation time tjj. In section VI we give final results for rn and in section VII we analyze the 
incorporation of the dilaton and calculate the expectation value of trF 2 . We close the paper with a discussion. 
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II. BOOST INVARIANT KINEMATICS 



An interesting kinematical regime of the expanding plasma is the so-called central rapidity region. There, as was 
suggested by Bjorken (l7j . one assumes that the system is invariant under longitudinal boosts. This assumption is in 
fact commonly used in realistic hydrodynamic simulations of QGP [2]. If in addition we assume no dependence on 
transverse coordinates (a limit of infinitely large nuclei) the dynamics simplifies enormously. 
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In order to study boost-invariant plasma configurations it is convenient to pass from Minkowski coordinates 
(i , x , x±) to proper-time/spacetime rapidity ones (r, y,x±) through 

.T°=rcoshy x 1 =rsinhy (I) 

The object of this work is to describe the spacetime dependence of the energy-momentum tensor of a boost-invariant 
plasma in J\f = 4 SYM theory at strong coupling. The symmetries of the problem reduce the number of independent 
components of to three. Energy-momentum conservation d^T^ = and tracelessness Tff = allows to express 
all components in terms of just a single function - the energy density s(t) in the local rest frame. Explicitly we have 
i 

T TT = e( T ) (2) 
T yy = -r 2 (V) + ± E (r)j (3) 

T xx = e(r) + \r^e{r) (4) 

Gauge theory dynamics should now pick out a definite function e(r). The aim of this paper is to find its behaviour at 
large proper-times, up to subsubasymptotic terms, and to interpret this behaviour in terms of parameters of second 
order hydrodynamics which we will describe in the next section. 



III. SECOND ORDER VISCOUS HYDRODYNAMICS 



The object of a hydrodynamic model is to determine the spacetime dependence of the energy-momentum tensor 
for an expanding (plasma) system. The simplest dynamical assumption is that of a perfect fluid. This amounts to 
assuming that the energy momentum has the form 

T^ u = (e + p)u^u v + pt]^ (5) 

where u M is the local 4- velocity of the fluid (u 2 = —I), e is the energy density and p is the pressure. In the case of 
TV = 4 SYM theory that we consider here Tff = and hence e = 3p. The equation of motion that one obtains from 
energy conservation in the boost-invariant setup is 

« £ + P 4e ra\ 
d T e = = --- (6) 

T 6 T 



whose solution is the celebrated Bjorken result 



T ■ ', 



(7) 



Once one wants to include dissipative effects coming from shear viscosity, the description becomes more complex. In 
a first approximation one adds to the perfect fluid tensor a dissipative contribution ?7(V Al u ly + V^u^), where rj is the 
shear viscosity of the fluid. The resulting equations of motion get modified to 



= + £ (8) 



Note that in the above equation the shear viscosity is generically temperature dependent cx T 3 in the Af = 4 case) 
and hence r dependent. In order to have a closed system of equations we have to incorporate this dependence through 

77 = A ■ (9) 

with A being some numerical coefficient. 

However this so-called first order formalism suffers from a number of problems. Firstly it is inconsistent with 
relativistic invariance (causality) - excitations may propagate at speeds faster than light. Secondly these equations 
suffer from some unphysical behaviour (see e.g. [161]). 

In order to cure these problems Israel and Stewart introduced a second order theory [l5| with an additional 
parameter - the so-called relaxation time Tn which can overcome the problems with causality. This theory has found 



applications in modelling heavy-ion collisions [18j |. The corresponding equations of motion, again in the Bjorken 
regime which we are considering here, are now 

4 e $ 

d T e = --- + - (10) 
At t 

r n 9 r * = -$ + ^ (11) 

where $ is related to the dissipative part of the energy momentum tensor ILj„ through $ = — T 2 H yy . Let us note that 
when Tn — > the equations reduce to the first order formalism case. In the above equations one has to keep in mind 
that the shear viscosity r\ is again given by @. Hence the only remaining independent parameter is therefore Tn- 

In the derivation of the second order viscous hydrodynamics from Boltzmann equations one finds that the ratio of 
Tn to j] is given in terms of the pressure p [l6[ 

Boltzmann = ^ = ^7 ^ 

2p 

This is the value usually used in viscous hydrodynamic simulations. Tha aim of this work is to determine the relaxation 
time in the strong coupling regime using the AdS/CFT correspondence. To this end let us parametrize the relaxation 
time at strong coupling as 

_ ^Boltzmann /-t o\ 

Tn = r ■ T n (13) 

where T B° ltzmann [ s defined by the expression (\V2\i and r is a numerical coefficient. The aim of this paper is to 
determine r. 



IV. ADS/CFT DESCRIPTION OF AN EXPANDING BOOST-INVARIANT PLASMA 

The procedure adopted in [8| to describe, using AdS/CFT, an expanding system of plasma in J\f = 4 SYM is (i) 
consider a family of possible behaviours of the spacetime expectation value of the energy-momentum tensor (T^„), 
(ii) for each of those (T M1/ )'s find the dual geometry which generically will be singular and (iii) use the criterion of 
nonsingularity of the constructed geometry to pick out the physical spacetime profile of {T^ v ) . 

For the plasma configuration considered in this paper, as described in section II, all possible possible boost-invariant 
profiles of (Tpv) can be expressed in terms of a single function e(r) which is just the energy density in the local rest 
frame. 

The construction of a dual geometry then proceeds as follows [19| . First we adopt the Fefferman-Graham coordinates 
(20| for the 5-dimensional metric 

2 ~g uv dx^dx v + dz 2 
ds = — (14) 

z z 

where the z coordinate is the 'fifth' coordinate while [i is a 4D index. g M „ is here a function both of the 4D spacetime 
coordinates and of the 'fifth' coordinate z. Then we have to solve Einstein equations with negative cosmological 
constant: 

E a(3 = R a/3 - -g a pR - 6g a (3 = (15) 
with a boundary condition for g M „ around z = 0: 

~9n V = + z A ~g { $ + ... (16) 
where the fourth order term is related to the expectation value of the energy-momentum tensor through 

N 2 

(2W = ^ (^) 



It is convenient to ignore the factor N 2 /(2ir 2 ) throughout the calculation and only reinstate it in the final result when 
e.g. writing the energy density or shear viscosity directly in terms of the temperature. 



For the boost-invariant geometries relevant here, the most general metric in the Fefferman-Graham coordinates 
takes the form 

ds 2 = -i (-e a ^dr 2 + e b ^T 2 dy 2 + e c ^dx 2 ± ) + ^ (18) 
z z \ J z z 

with three coefficient functions a(z,r), b(z, r) and c(z,t). The energy density defining the physics is simply related 
to the boundary asymptotics of the a(z,r) coefficient: 

eT )=-Iim-^ (19) 

Of course it is too difficult to explicitly perform this construction for arbitrary functions e(r). What one does in 
practice is to perform an expansion of e(r) for (large) proper-times r and determine one by one the subsequent terms 
in the expansion. In the following section we will review the results obtained so far in carrying out this program 



REVIEW OF THE DUAL GEOMETRY UP TO O r 



In [| it was shown that in order to study large proper-time limit of the metric one is led to introduce a scaling 
variable 

v = 4 (20) 

T3 

and take the limit r — > oo with v fixed. In order to study subleading terms in the metric we have to perform an 
expansion around this limt. 

Let us now expand the metric coefficient functions a(z,r), b(z,r) and c{z,t) in a series of the following form: 

a(z,r) = a (v) + a x {v)\- + a 2 {v)\ + a 3 (v)\ + ... (21) 

73 T3 T z 

and similar expressions for the other coefficients. The motivation for choosing such specific powers of r is twofold. 
Firstly, on the gauge theory side once we have viscous hydrodynamics with r\ oc T 3 , the asym pto tic expansion of 
the energy density e[r) exactly corresponds to the above decomposition of the metric coefficients j2a] (|2T|) . Secondly, 
one can show directly from the gravity side without assuming anything on viscous dynamics in gauge theory, that a 
correction at order l/ra has to occur. Namely suppose that we start from the leading order solution with a generic 
first correction 

a{z,T) = a(v)+a r (v)— + ... (22) 

T 

Then we find that the square of the Riemann tensor is nonsingular at the leading and l/r r orders but will always 
have a singularity at order I/ts . There will also be a singularity at order l/r 2r . The only possibility of obtaining a 
nonsingular geometry is that 2r = | and that these two singularities cancel, which is indeed what occurs. This fixes 

the power of the first subleading correction to be l/ri which is exactly what is expected for the form of corrections 
due to shear viscosity with the coefficient behaving like n oc T 3 . 

The procedure is now to insert the expansion (|21[) into Einstein's equations and solve them order by order. At 
each order a new integration constant (free parameter) will occur and we will determine it by requiring that a similar 
expansion of the square of the Riemann tensor will be nonsingular i.e. 

m 2 = R^ af3 R^ a0 = R (v) + Rx{v)\ + R 2 {v)\ + R 3 (v)^ + ... (23) 

with all Ri(v) being nonsingular. 

In order to present the solutions it turns out to be convenient to define 

6 j («) = 6 i (») + 2c i («) (24) 



The leading order solution found in [8j is 

(1 - ^ 4 /3) 2 
5o(«) = log(l + « 4 /3) 3 

co(«) = log(l + « 4 /3) (25) 
and the resulting Dl 2 coefBcient is nonsingular: 

8(5v 16 + 60v 12 + 1566?/ + 540u 4 + 405) 

Ro = (26) 

From this expression, using the similarity with the static black hole metric we may read off the temperature from the 
position of the horizon to obtain at this order: 

T(t) = -^K (27) 

0*7TT 3 

We will use this expression later in the paper. The first subleading correction was found in [13] and reads in our 
conventions 

(9 + v 4 )v 4 

a M = 2r io g _ yS 

h(v) = -6^o^ 

<*(«) = ^o^-r/olog^ (28) 

where 770 is an undetermined integration constant (which has the physical interpretation as the coefficient of shear 
viscosity) . However it turns out that r/o is undetermined at this order from nonsingularity of d\ 2 since 



41472(v 4 - 3)v 8 



and is nonsingular for any value of t/q. Therefore in order to fix 770 one has to go one order higher. This was done in 
[§] with the result 

(9 + 5v 4 )v 2 ^(9 + v 4 )v 4 ,(-1053 - 171w 4 + 9u 8 + 7w 12 )f 4 1 , VS - v 2 3 2l 3 - v 4 



6(9-u 8 ) 36(9-« 8 ) ' u 3(9-tj 8 ) 2 4\/3 V^ + v 2 2 ,u b 3 + v 4 

v 2 v 4 2 (39 + 7v 4 )v 4 1 VS~v 2 3 2 3 -w 4 



2(3 + « 4 ) 12(3 + « 4 ) ,u (3 + « 4 ) 2 4^3 ^^3 + ^ 2 ,u b 3 + 7/ 
u 2 (9 + w 4 ) u 4 2 (-9 + 54?/ + 7i; f > 4 1 \/3- 



2 



C2{V) 144V3 + 6(9- « 8 ) +C 36(3 + ^) % " 3(3 + W 4 )(9-^ 8 ) + 4^3 l ° g VI + v 2 + 



1 (/~t , 2\ i 3-« 4 1 / V3-v\ (V3-v 2 )(V3 + v 2 r .. ( (V3-^ 2 ) 2 ' 
" 36 + 66 % ) log 3^ + ^ (log log h 2 I "Tyf^ 



Here C is a new free integration constant appearing at this order and li 2 is the dilogarithm function. A calculation of 
the V\ 2 coefficient at this order gives 



576(v 4 - 3)u 8 6912(5w 24 - 60tj 20 + 2313tj 16 - 6912w 12 + 26487u 8 - 18468w 4 + 13851)*; s rj 2 
2 ~ (3 + tj 4 )5 + (3-i; 4 ) 4 (3 + ?j 4 ) 6 

4608(5v 16 + 6v 12 + 162v s + 54w 4 + 405)u 10 
(3-v 4 ) 4 (3 + v 4 ) 5 



(31) 



We see that there is a fourth order pole singularity. It may be cancelled when the viscosity coefficient 770 takes the 
value 



which reproduces, in the boost-invariant expanding setup, the exact viscosity coefficient of N = 4 SYM calculated in 
the static case in 

We see however, that at this order the new integration constant C is still undetermined. By analogy with the case 
of r]o we may expect that it will be fixed at the next order of the series expansion. Before we proceed to do this let 
us first discuss the physical interpretation of this coefficient. 



by 



The physical interpretation of C 

The energy density extracted from the metric expanded in the scaling limit up to O y~^J through (|f 9p is given 



Let us now suppose that this behaviour can be described by the equations of second order viscous hydrodynamics 
(fTClf -Qf f jl with some parameters r\ and tjj- We will use this hypothesis to extract these parameters [2 7] which may 
then be used in hydrodynamic simulations for generic 3+1 evolving plasma systems. 

For ease of computation we note that equations (fl"U|) - (fTTj) can be divided out by the factor N 2 /(2ir 2 ), thus we can 
drop this factor from e, rj and <£>. Let us now insert $ from (fTT)|) into (fTTj) and obtain a differential equation expressed 
sorely in terms of e(r) and the numerical coefficients A and r. 

We may now extract these numerical coefficients by plugging in the energy density from the AdS/CFT calculation 

into the resulting equation. 
The result is 

A = vo (34) 
-C-66772 



324^ 



(35) 



A is precisely equal to r/o which is not unexpected since in the leading order e(r) = 1/t3 and hence r\ = t]q/t may be 
written as r/ — Tfo£%. The value of 770 is determined from nonsingularity to this order [j| to be given by (|32|) . The new 
information is r defining the relaxation time rn- Plugging in the value (|32p for 770 we may express r in terms of C: 

r= -i-i^i (36) 

Hence the relaxation time rn will be found once we know C (we note that thus we must have C < 0). Again, as was 
the case with 770 one has to go one order higher in order to determine the value of the coefficient. We will proceed to 
do it in the following section. 



VI. DETERMINATION OF C AND THE RELAXATION TIME m 

It is straightforward to obtain equations for the third order metric coefficients as (v), b 3 (v) and Ca(v) using a 
computer algebra system, albeit these equations appear to be prohibitively complex at hrst sight. 

However one can notice that the zz component of the Einstein equations E zz — (|15p gives an equation just for 
the combination 

d' 3 (z) ee a' 3 (v) + b' 3 (v) ee a' 3 (v) + ^(t,) + 2c' 3 (v) (37) 

where the prime stands for the ordinary derivative. Then expressing b' 3 (v) in terms of the other functions and plugging 
the result into the TT component E TT = of Einstein equation yields an equation purely for the derivative a' 3 (v). So 
at this stage we have analytic expressions for a' 3 (v) and (say) b' 3 (v). 

Fortunately, the expression for the O (t~ 2 ) coefficient R 3 of 9\ 2 can be expressed sorely in terms of a' 3 (v), b 3 (v) 
and their derivatives. The expression is quite lengthy so we do not give it here [2l|. Substituting a' 3 (v) and b 3 (v) into 
i?3, we find again that there is a fourth order pole in R 3 at v = 3* which is canceled exactly when 

C= 6l0g2 - 17 (38) 
y/3 



For this value the poles of lower orders are also canceled. However a new feature arises here, namely there is a leftover 
logarithmic singularity: 



i?3ii r _ 6io 8 2-i7 = finite + 8 • 2^3* log(3 - w 4 ) (39) 

We will show in the next section that if we turn on the dilaton this singularity may be canceled without modifying 
the value of C determined above. Let us now determine the coefficient of the relaxation time r from (|36|) . We find 
that 

r = 1(1 -log 2) ~ 0.034 ~ i_ (40) 

which shows that the relaxation time of Af = 4 SYM at strong coupling is much smaller than one would expect. 
In particular the second order viscous hydrodynamics seem to be much closer to the first order behaviour. If we 
express the relaxation time in terms of the proper time and then relate it to the temperature through the leading 
order expression (|27[) . we obtain 

l-log2 

Tn = "W (41) 



VII. THE DILATON AND trF 2 



In this section we will show how to cancel the remaining logarithmic divergence in i?3 by turning on the dilaton 
field. On the gauge theory side, this means that we generate a nonzero expectation value for tr F 2 i.e. electric and 
magnetic modes are no longer equilibrated. Physically this may well happen since we are considering a regime where 
dissipative effects are important. 

If we include the dilaton </>, the equation of motion become those of a coupled Einstcin-dilaton system which read 
(in the Einstein frame): 



1 

2 l 

dnWggTW) = (43) 



R^ + tg^u = ^<j>d v <p (42) 



In analogy to the metric coefficients (l2"Tj) we would like to make a large r expansion of the dilaton. This requires a 
judicious choice of the powers of r appearing in that expansion. We fix those powers by the following argument. 
Suppose that 4>(z, t) ~ <f)(v)T~ r . Then one can check using (|4"2"|l that the dilaton source term will contribute to the 

The above considerations suggest the following expansion of the dilaton field: 



z ,r) = Y J Uv)\ (44) 



The i = component is absent as has already been checked in [l2| ■ 

If a spacetime has a nontrivial dilaton profile, in string theory one has two distinguished metrics - the Einstein 
frame metric considered above, and the string frame metric defined as 

9f v ma = e* V (45) 

which is the natural metric from the point of view of strings propagating in the spacetime. Therefore a natural 
question arises to which of the two metrics we should apply the nonsingularity criterion. A natural guess would be 
that %*?{g 8 av %n9 ) should be nonsingular, we will however for the moment keep both possibilities open. 

Let us now go back to the expansion (|4"4")) . The only terms that would contribute to the metric coefficients [28] 
a\(v), a,2(v) and ci3(v) are cj>i(v), 4>2( v ) an( i <p3 ( v )- So only these terms would influence 9\ 2 in the Einstein frame to 
the order considered in the previous sections. If on the other hand we would require the nonsingularity of V\ 2 in the 
string frame, then we must consider in addition <f>4(v), <f>s(v) and (f>e (v). 

Finally let us give one gauge theoretical interpretation of a nontrivial dilaton profile. It might give rise to an 
expectation value for tr F 2 through [22|, [23| 



JYM 27r2 Z ^ Q 



In the following we will consider various possible leading behaviours of the dilaton according to the expansion 
and consider 5H 2 in string frame (and in Einstein frame where relevant) 

at O (t -1 / 3 ). At this order the leading behaviour of can be easily obtained as 

1 („)=Mogf^J ( 4? ) 

3 + V 4 

In the string frame SRj tri will have a new contribution at order r -1 / 3 . This contibution has a piece which is 
proportional to fclog(3 — v 4 ). Therefore we have to set k = and hence 4>i(v) = 0. Let us note that the natural 
appearance of a logarithmic divergence gives a hope that indeed the dilaton might cancel the leftover divergence in 
(155)) . We will see below that this indeed happens. 

In the Einstein frame the calculations are more complicated. 0i will modify the metric coefficients a\{v), b±(v) and 
c\(v). Then the Einstein frame fH 2 will also get modified and (f2"9"f will be replaced by 

41472(v 4 - 3)i> 8 1152(13w 12 - 99i; 8 + 27v 4 + 27)v % 2 

Rl= (3 + ^)5 (3 + w 4 )5( 3_ w 4 )2 k ( 48 ) 

Requiring nonsingularity again gives k = and hence 4>i{v) = 0. 

at O (t~ 2 / 3 ). The leading part 02 (v) has a functional form identical to (|47|) . Again the string frame 9i 2 trin9 will 
have a logarithmic singularity but now at order r~ 2 / 3 , hence 02 {v) has to vanish at this order if we assume string 
frame nonsingularity. 

In the Einstein frame, the coefficients d2{v), 62 (v) and C2(u) will get modified. EH 2 will still have a fourth order pole 
at O (t ~ 4 / 3 ) which will be canceled by taking (|32p but an additional second order pole proportional to k will persist. 
Hence also requiring nonsingularity in the Einstein frame will lead to 02 (v) = 0. 

at O (r -1 ). Here again in the string frame one has a logarithmic singularity hence <f>s(v) = 0. The calculations in 
Einstein frame are very tedious but also rule out a nonvanishing 03 (v). 

at O (t~ 4 / 3 ) . At this order, the metric coefficients are not modified, hence it makes sense to consider only the 
string frame ^1 tr i ng - A logarithmic singularity which appears at order r -4 / 3 requires 4>a{v) = 0. 

at G(r- 5 /3). A gain a logarithmic singularity requires 05 (v) = 0. 

at O (t -2 ). This is the relevant order at which a cancelation of the logarithmic singularity may occur. At this 
order 06 has again the form 

6 ( V ) = fclog|-4 (49) 
3 + w 4 

The string frame EH 2 triTlg at O (t~ 2 ) will get modified by 

pSt ™ 3 D , 8(5^ 16 + mv 12 + 1566z; 8 + 540^ 4 + 405) 3 - ^ 4 

i? 3 = i? 3 - fc log ——j (50) 



where R3 is the Einstein frame coefficient calculated earlier. Therefore we see that by a suitable choice of k we may 
exactly cancel the leftover logarithmic singularity without changing the value of C obtained above which came from 
canceling a fourth order pole. Indeed performing the expansion of R s ^ rm9 at u = 3 1 / 4 we get 



L 3 

which gives 



R string = + g , tftf _ mfc ^3 _ ^ ^ 



0aW^log|^ (52) 



Now we may evaluate the resulting expectation value of (trF 2 ) from (|4U|) : 



1 ^ (53 , 



4<? 2 M N ' 2tt 2 7-3^ 
We see that this expectation value is negative which signifies that electric modes become dominant 



VIII. DISCUSSION 



In this paper we have considered the subasymptotic proper-time evolution of an infinite expanding boost-invariant 
plasma in M = 4 SYM at strong coupling. We have used the AdS/CFT correspondence to construct dual geometries 
and used the criterion that the unique nonsingular geometry corresponds to the physical evolution of the plasma. 
Then the proper-time dependence of the energy density e(r) can be read off from the form of the 5-dimensional 
metric. 

We have found subleading terms of e(r) which are sensitive to the relaxation time rn appearing, in addition to the 
shear viscosity 77, as a new element in second order viscous hydrodynamics. This served to determine the relaxation 
time pUlt - tETTj) which turned out to be about thirty times shorter than the one expected from Boltzmann kinetic theory 
estimates (|12[) . Thus second order hydrodynamics appears to be much closer to the ordinary first order formalism. 

In addition we have found that canceling a remaining logarithmic divergence requires turning on the dilaton field 29]. 
A nonvanishing dilaton field in turn implies a nonzero expectation value (|53[) for tr F 2 with the proper-time scaling 
T -io/3 This means that electric and magnetic modes are not exactly equilibrated. 

It would be very interesting to understand from a microscopic perspective the very short relaxation time. In 
particular it would be interesting to understand the weak coupling corrections to the Boltzmann value T ^ oltzmam \ 
Another intriguing question would be to estimate by other means the r dependence of (tr F 2 } . On the strong coupling 
side it would be interesting to understand the dual metric from the point of view of dynamical horizons in general 
relativity [2a ] and analyze the relevant thermodynamics. 
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